Scattering of two-dimensional solitons in dipolar Bose-Einstein condensates 
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We analyze the scattering of bright solitons in dipolar Bose-Einstein condensates placed in uncon- 
nected layers. Whereas for short-range interactions unconnected layers are independent, a remark- 
able consequence of the dipole interaction is the appearance of novel nonlocal interlayer effects. In 
particular, we show that the interlayer interaction leads to an effective molecular potential between 
disconnected solitons, inducing a complex scattering physics between them, which includes inelastic 
fusion into soliton-molecules, and strong symmetric and asymmetric inelastic resonances. In ad- 
dition, a fundamentally new 2D scattering scenario in matter-wave solitons is possible, in which 
inelastic spiraling occurs, resembling phenomena in photorrefractive materials. Finally, we consider 
the scattering of unconnected ID solitons and discuss the feasibility in current on going experiments. 
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Up to very recently, typical experiments on ultra cold 
gases involved particles interacting dominantly via a 
short-range isotropic potential, which, due to the very 
low energies involved, is fully determined by the corre- 
sponding s-wave scattering length. However, recent ex- 
periments on cold molecules 0, Rydberg atoms Q, and 
atoms with large magnetic moment 0, open a faTcinat- 
ing new research area, namely that of dipolar gases, for 
which the dipole-dipole interaction (DDI) plays a signif- 
icant or even dominant role. The DDI is long-range and 
anisotropic (partially attractive) , and leads to fundamen- 
tally new physics in condensates 0, 0, @ , degenerated 
Fermi gases and strongly-correlated atomic systems 
@. It leads to the Einstein-de Haas effect in spinor con- 
densates and may be employed for quantum com- 
putation lfj, and ultra cold chemistry Recently, 
time-of-fiight experiments in Chromium have allowed for 
the first observation ever of dipolar effects in quantum 
gases |T^ . 

Interestingly, the physics of short-range interacting 
Bose-Einstein condensates (BECs) at low temperatures is 
given by a nonlinear Schrodinger equation (NLSE) with 
cubic local nonlinearity, similar to the one appearing in 
other systems, in particular in Kerr media in nonlinear 
tics. In ID, nonlinearity allows for solitonic solutions 
which have been observed in BEC |LJ|. However, 
bright solitons are unstable in 2D and 3D. In periodic 
potentials multidimensional discrete solitons are possible 
|l5| . but they do not move in a multidimensional way, 
although the use of optical lattices has been proposed to 
move 2D and 3D discrete solitons along a free direction 
[lif . Other interesting possibility to stabilize high dimen- 
sional solitons is to use Feshbach resonances to manage 
spatially and/or temporally the scattering length [l7j . 

Due to the DDI, a dipolar BEC is described by a NLSE 
with nonlocal cubic nonlinearity [J opening an in- 
teresting cross-disciplinary link between BEC and other 
nonlocal nonlinear media, as e.g. plasmas |l8j . where 
the nonlocal response is induced by heating and ioniza- 
tion, and nematic liquid crystals, where it is the result of 



long-range molecular interactions [19| . Nonlocality plays 
a crucial role in the physics of solitons and modulation in- 
stability [U mil. 

In particular, any symmetric nonlo- 
cal nonlinear response with positive definite Fourier spec- 
trum has been mathematically shown to arrest collapse 
in arbitrary dimensions |2lj . Indeed, multidimensional 
solitons have been experimentally observed in nematic 
liquid crystals . Recently we showed that under real- 
istic conditions, 2D solitons may be generated in dipolar 
BEC [23^ . However, the anisotropic character of the DDI 
violates the conditions of Ref . [2l| , and as a consequence 
a stability window occurs, rather than a stability thresh- 
old for a sufficiently large dipole strength. In Ref. [2^] 
we briefly studied the scattering of 2D dipolar solitons, 
which is inelastic j^, contrary to the ID solitons in lo- 
cal NLSE, due to the lack of integrability |25|. However, 
the analysis of the inelastic scattering was largely com- 
plicated by the spatial overlapping of the solitons. 
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Figure 1: Schematic representation of the system considered 
in this Letter. 

In the following, we consider dipolar BECs placed in 
disconnected layers (Fig. For the case of short-range 
interactions, the two layers are independent if the hop- 
ping is suppressed. A remarkable consequence of the DDI 
is that unconnected layers become coupled due to non- 
local density-density interactions, leading to interesting 
interlayer effects, as e.g. the possibility of a BEC of fila- 
ments, recently discussed in Ref. |2g. In this Letter, we 
analyze the rich physics introduced by interlayer effects in 
the nonlinear properties of dipolar BECs. In particular, 
we show that this interlayer interaction leads to an effec- 
tive molecular potential between fully disconnected soli- 
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tons, allowing for a complex scattering physics between 
them. This novel physics includes inelastic fusion into 
excited soliton- molecules for sufficiently slow solitons, as 
well as strong symmetric and asymmetric inelastic reso- 
nances for intermediate velocities. In addition, we discuss 
a fundamentally new 2D-scattering scenario in matter- 
waves, showing that inelastic soliton spiraling similar to 
that observed in photorrefractive materials |22j,l22| is pos- 
sible in dipolar BEC. Finally we consider the scattering of 
ID dipolar solitons in unconnected wires, and comment 
about observability in on-going experiments. 

In the following, we consider a dipolar BEC transver- 
sally confined in the z-direction by a two-well potential, 
with wells located at z = ±zo, and separated by a suf- 
ficiently large potential barrier that prevents tunneling 
between them. At each well the ^-confinement is ap- 
proximated by an harmonic potential of frequency uj z , 
whereas there is no confinement on the xy-plane. We 
consider in each well a BEC of N particles (a more gen- 
eral case will be analyzed later on) with electric dipole 
d (the results are equally valid for magnetic dipoles) ori- 
ented in the z-direction by a sufficiently large external 
field, and that hence interact via a dipole-dipole poten- 
tial: Vdif) = gd{l -3 cos 2 6>)/r 3 , where g d = aNd 2 /4neo, 
with eo the vacuum permittivity, 9 the angle formed by 
the vector joining the interacting particles and the dipole 
direction, and — 1/2 < a < 1 a tunable parameter by 
means of rotating orienting fields [2^| . At sufficiently low 
temperatures our system is described by the following 
two coupled NLSE with nonlocal nonlinearity 0: 
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where j = ±1 is the layer-index, ^ j are the wavefunctions 
at each well, Uj(z) = muj 2 (z + jz ) 2 /2, J |%(f, t)\ 2 df = 
1, and g = Aith 2 aN / m characterizes the contact interac- 
tion, with a the s-wavc scattering length. In the following 
we consider a > 0, i.e. repulsive short-range interactions. 

We assume a 2D dynamics in each well. This approxi- 
mation demands that the corresponding chemical poten- 
tial \x <C fko z . In that case, ^j(f) ~ ipj{p)ipj(z) with 
¥>j{z) the ground-state wave-function of the harmonic 
oscillator in the layer j. Employing this factorization, 
the convolution theorem, the Fourier transform of the 
dipole-dipole potential, Vd(k) — (4ir/3)(3k 2 /k 2 — 1), and 
integrating over the z-direction, we arrive at a system of 
two coupled 2D NLSE: 
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(h 3 (k p )F {k p l z ,0) + h-j{k p )F (k p l Zl 2z /l z )^ (2) 
where I 2 = h/muj z is the harmonic-oscillator length, hj 



is the Fourier transform of \ipj(p)\ 2 arj d F(\/2k, V2X) = 
2e- A ' - (3 V^fce fc2 /2) [e- 2kx erfc(fc - A) + e 2feA erfc(fc + A)] , 
with erfc(a;) the complementary error function. Using 
Eqs. J5J), we study numerically the equilibrium properties 
and the dynamics of the unconnected 2D solitons. 

In Ref. jl^l we showed that if gd ^ 0, a stable soliton 
may appear if (3g < 3(1 + g/2ir)/2 < —2/3g, where g = 
g j 1 \f2ihio z l\ 1 and (3 = gd/g- Hence [3 must be sufficiently 
large and negative , which demands tunability (a < 0). 
For Na/l z ^> 1, stable solutions appear for \f3\ > 3/87T ~ 
0.12. As discussed in Ref. [23, the DDI may destabilize 
the solitons if they become 3D. In the following, we shall 
restrict ourselves to the 2D regime |23j . 

We introduce a variational formalism which allows us 
to study the equilibrium properties and dynamics of the 
two solitons. We consider a Gaussian Ansatz 301: 
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where A is the normalization factor, {xq, yg} is the soliton 
center, w XtV the soliton widths, and from the continuity 
equation we obtain a v = mrjo/h, (3^ = mw ri /2huj ri . The 
variables xo,yo,w Xl w y are time-dependent. The center 
of mass motion is an independent degree of freedom and 
it can be decoupled. Without loss of generality, it has 
not been included in the variational Ansatz. Introducing 
© into the corresponding Lagrangian |22| , we obtain the 
following set of equations of motion 
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where q^y = xo,yo,w x ,w y are the dynamical variables, 
and 
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+ V, (5) 



is the potential energy, that includes the dipolar interac- 
tion term 

V = J±- I d kU%~^e-^^-^l 2 - k >V 2 
\2tt z J \ k l J 

(1 + cos(2fc;;Xo) cos(2k y yo) cos(2fc z zo)), (6) 

that couples the unconnected solitons. Hence, the prob- 
lem reduces to an effective particle in the potential U. 
Since we have set gd < 0, the soliton-soliton poten- 
tial is maximally repulsive for solitons on top of each 
other, becoming attractive at a given distance. A soliton 
molecule is thus possible at the minimum of U, which 
we have found by means of a Powell-minimization pro- 
cedure, obtaining results that compare well with our 
imaginary-time simulations of Eqs. (0. For point- like 
solitons the soliton-soliton potential would be V po i n t oc 
— (xq — 2zq)/(x 2 + Zq) 5 / 2 , with a minimum at xq = 2zq, 
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Figure 2: Numerical (crosses) and variational (solid) results 
for Afc/fco (Afe = ko — k(t — > oo)) as a function of the initial 
momentum fco4> for 2o = 34, p = 200, /3 = —0.2. Inset: 
Numerical results with zq = 41 z (solid) and zo = 5Z Z (dotted). 
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Figure 4: Numerical result for the soliton trajectory during 
spiraling fusion in a 2D soliton scattering, for the case g = 
200, /3 = -0.2, 2 = 34, fc 4 = 0.01a, x = 304, y = 104 . 



Figure 3: Variational results for the evolution of the averaged 
width w p = sj w% + w% (normalized to its initial value) for 
soliton 1 with Ni = 0.6iV (dotted) and 2 with N 2 = IAN 
(solid), with g = 200, f3 = -0.2, z = 34, and fe 4 = 0.03. 
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Figure 5: Width of a ID soliton for a;o = 34, g/2itfiu! p l p = 
25, /3 = 0.28, and k l p = 0.05 (solid), 0.20 (dashed), 0.35 
(dotted-dashed). uo p (l p ) is the transversal oscillator frequency 
(length). The time has been re-scaled for comparison. 



a value much smaller than that obtained from our vari- 
ational or numerical calculations, showing the relevance 
of the spatial extension of the solitons. 

In the following, we consider the scattering of solitons. 
We first discuss the case where the relative velocity is par- 
allel to the vector connecting the centers of mass of the 
two solitons (yo = 0). We have studied the scattering 
for different initial soliton velocities both by direct nu- 
merical simulations of Eqs. (J2J) and by determining the 
evolution of {xo,w x ,w y } in our variational calculation. 
Fig- shows the relative variation of the soliton momen- 
tum as a function of the initial momentum. As expected, 
for sufficiently large initial velocities the scattering may 
be considered as elastic. For sufficiently low velocities 
the initial kinetic energy of the solitons is fully trans- 
formed during the inelastic scattering into internal soli- 
ton energy, and the initially independent solitons become 
bounded into an excited molecular state (fusion). 

Interestingly, the inelastic losses do not increase mono- 
tonically for decreasing velocities, but on the contrary 
show a pronounced resonant peak at intermediate veloc- 
ities (Fig. Eland its inset). This effect is motivated by a 
resonant coupling to internal soliton modes, which leads 
to a dramatic enhancement of the soliton widths after 
the collision, and eventually to the destruction of the 



solitons. We stress that this is only possible because in- 
ternal modes of the 2D soliton are at rather low energies, 
well within the inelastic regime. If the interlayer distance 
is increased, the inelastic losses are as expected reduced, 
but an even more complicated structure of resonances is 
then resolved (Fig. inset). 

The previous analysis can be extended to asymmetric 
configurations, where the solitons have different number 
of particles. In this case, asymmetric inelastic scattering 
occurs, as shown in Fig. [21 where a resonance appears 
just for the solitons having the smaller binding energy. 

Interestingly, the possibility of generating stable 2D 
solitons in dipolar gases allows for a completely new sce- 
nario for soliton scattering in cold gases, namely a truly 
2D scattering, in which the scattered solitons present a 
relative angular momentum around the scattering cen- 
ter (Fig. 01 . Similar scattering regimes as those dis- 
cussed above are also possible in the 2D scattering. How- 
ever, the relative angular momentum during the collision 
leads, for the 2D scattering case, to a spiraling motion 
for sufficiently low incoming velocities, as a consequence 
of the inelastic fusion of the solitons, which stabilize in 
a rosetta-like orbit around each other (Fig. 0. The 2D 
spiraling links the physics of dipolar BECs to that of pho- 
torrefractive materials, where soliton spiraling has been 
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proposed (2^ and experimentally observed |2s| . 

Finally, we consider ID BECs placed at neighboring 
ID-sites (x = ±xq) of a 2D lattice, with the dipole ori- 
ented along the site axis. Following similar scaling ar- 
guments as above, it is possible to show that for g > 0, 
a bright soliton is possible if p > 3/4-7T. Although, of 
course, some 2D features are missed in ID, it is indeed 
possible to observe inelastic processes also in ID solitons. 
This is particularly relevant for current experiments in 
52 Cr, since no tuning is necessary for ID solitons, eas- 
ing very significantly the experimental requirements. For 
52 Cr a Feshbach resonance is necessary to satisfy the pre- 
vious condition, but Feshbach resonances are well char- 
acterized and accessible [3lj . In Fig. we depict an ex- 
ample of the dynamics of the soliton width in and out 
of resonance, which clearly shows a resonant (although 
non-destructive) behavior for intermediate velocities. 

Summarizing, interlayer effects are a fundamentally 
new feature introduced by the DDI in dipolar gases 
placed in unconnected layers of an optical lattice. These 
effects may have remarkable consequences, as e.g. the 
formation of a BEC of filaments 26]. In this Letter, we 
discussed the rich physics introduced by interlayer effects 
in the nonlinear properties of dipolar BECs, and in par- 
ticular in the scattering of unconnected solitons. The 
DDI induces an inelastic soliton-soliton scattering, that 
for low relative velocities, leads to the inelastic fusion 
into a soliton molecule. Interestingly, the inelastic losses 
do not increase monotonically for decreasing relative ve- 
locities, but on the contrary show strong resonances at 
intermediate velocities, at which, after interacting, the 
soliton widths are strongly modified, eventually leading 
to soliton destruction. This effect appears, because, due 
to the relatively low excitation frequencies of the solitons, 
a resonant coupling between incoming kinetic energy and 
internal soliton modes is possible for low relative veloc- 
ities well within the inelastic regime. We have shown 
that a similar effect should be observable in ID geome- 
tries, where the experimental requirements may be easily 
fulfilled in on-going Chromium experiments. Finally, we 
have considered the 2D scattering of dipolar solitons, a 
unique possibility offered by the dipolar interactions in 
cold gases. We have shown that due to the combina- 
tion of inelastic trapping and initial angular momentum 
a spiraling motion is possible, offering fascinating links 
to similar physics in photorrefractive materials. 
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